Modelling unidirectional liquid spreading on slanted microposts 
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A lattice Boltzmann algorithm is used to simulate the slow spreading of drops on a surface patterned with 
slanted micro-posts. Gibb's pinning of the interface on the sides or top of the posts leads to unidirectional 
spreading over a wide range of contact angles and inclination angles of the posts. Regimes for spreading in 
no, one or two directions are identified, and shown to agree well with a two-dimensional theory proposed in 
Chu, Xiao and Wang fll. A more detailed numerical analysis of the contact line shapes allows us to understand 
deviations from the two dimensional model, and to identify the shapes of the pinned interfaces. 

I. INTRODUCTION 

Unusual wetting and spreading properties of drops on natural and artificial surfaces can often be explained by the micro- 
structure of the substrate Q. For example drops spreading on superhydrophobic surfaces patterned with micron-scale ridges 
reach an elongated final state fT/?], and a drop imbibing into a surface patterned with posts can form a faceted final configura- 
tion, which reflects the symmetry of both the lattice and of the posts themselves |5|. 

The physics behind this behaviour was first described by Gibbs, who pointed out that an interface can pin on the edge of a 
post over a range of angles, as illustrated in Fig. [T] The pinning occurs because there is a free energy penalty to the interface 
moving away from the edge in either direction as it would then have to form an angle with the adjacent surface which differs 
from the equilibrium contact angle. A pinning strength that depends on the lattice direction leads to the drop having one or more 
preferred directions of motion, and hence anisotropic drop movement and shapes. 

In this paper we focus on unidirectional drop motion: where the symmetry of the underlying surface structure can pick out 
one easy direction of spreading ISHT)- Such surfaces occur naturally, for example, the unidirectional motion of droplets on 
butterfly wings results from their ratchet-like structure f8^ '91, and rye-grass leaves shed water in a preferred direction ifTOl 
due to the asymmetric contact angle hysteresis. Microfabricated surfaces that lead to uni-directional motion are a very recent 
development. Unidirectional spreading has been observed on bent silicon micro pillars while other authors ifTTllTSl were 
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FIG. 1 : Gibbs' pinning on the comer of a post. The interface (blue line) remains pinned to the post over the range of angles indicated by as 
there is a free energy barrier to its moving in either direction. 
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able to transport droplets on vibrating ratchet structures. Similar results where obtained using PDMS replicas of the naturally 
occurring asymmetric micro-texture |10| taken from rye grass, and anisotropic hysteresis was observed on printed ratchetlike 
surfaces |fT3l . Despite the different material and geometries employed, the unidirectional liquid motion consistently reflects 
the asymmetry of the substrate on microscopic length scales. This highlights the importance of understanding the underlying 
physical phenomena involved. 

In this paper, we use a two-phase lattice Boltzmann algorithm to model imbibition on an hydrophilic surface patterned with 
slanting posts, varying the contact angle of the substrate and the tilt angle of the posts. For a range of contact angles we observe 
a single, preferred spreading direction as observed in recent experiments |1 1. The results are in good qualitative agreement with 
a two-dimensional model of the uniaxial spreading proposed in 1 1 ], and enable us to describe the corrections to the model for a 
three dimensional geometry. By visualizing the shape of the contact line we describe in detail the mechanisms through which 
the interface pinning and de-pinning occurs. 



II. THE MODEL LIQUID 
A. Governing equations 

To model a two phase system interacting with a surface, we apply a diffuse interface scheme. The thermodynamic state of 
the fluid is described at every point x and time t by its density p{x, t). The free energy of the system, ^, is taken as a Landau 
double-well potential with the addition of a derivative term representing the surface tension, and a surface contribution of the 
form proposed by Cahn l,14J : 
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The first term in the integrand of ([T]) is the bulk free energy density ifTSll 
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where , Pc , tw and /3 are, respectively, the critical density, critical pressure, reduced temperature and a free parameter 
controlling the density difference between phases and v — ^^p^ is a normalised density. This potential leads to two equilibrium 

bulk densities — pc (l ± y//3Tw)- Pb is a Lagrange multiplier constraining the total mass of fluid, while the third term in the 
free energy expression is an interface energy cost, tunable through the parameter k, associated with density gradients. It allows 
for solutions with a diffuse interface between phases, with surface tension 7 and width x- 
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The final term in Eq. ([T]l is the surface contribution to the free energy ^E*. When is minimised this gives the boundary condition 
d±p = —ps/i^ which fixes the value of the density at the solid surface. The Young contact angle at the surface is related to the 
surface chemical potential by ifTSl 



Ps = 2/3Tvi/\/2pcKsign ("^ ~ ^ 



cos |(l - COS I 



a — arccos(sin'^ 9). (4) 

In our simulations, the main parameters are set as follows: k = 0.01, Pc = 0.125, /3 = 1, tw = 0.3, pc = 3.5. The 
corresponding surface tension is 7 = 0.04 and the surface thickness is x = 0.9 (in simulation units). 

The hydrodynamics of the fluid is described by the continuity and Navier-Stokes equations: 

dtp + da{pua) ^Q, (5) 

dtipUa) + disipUaU/s) = 

-dpPaf} + dp {pi] [dpua + daup] + pXSa/sd^u^) , (6) 

where u is the fluid velocity field and 77 and A are the shear and bulk kinematic viscosities, respectively, dt represents a time 
derivative and da, dp spatial derivatives (Einstein summation convention is assumed). The connection between the thermody- 
namic and fluid dynamic of the system arises through the pressure tensor P, which is derived from the free energy ([T]i. Equations 
(|5]l and (j6]l are solved using the Lattice Boltzmann method. Details of the implementation can be found in [9l ll5( . 
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FIG. 2: (a) Diagrams to contrast the quasi-2D and full 3D geometries used in the simulations. The red and blue lines represent typical interface 
positions, (b) Side view of the posts showing the interface geometry corresponding to the threshold for movement which is assumed in deriving 



Eq. j7}. 



B. Geometry 



Each post has a square cross section of dimensions w, which was typically chosen to be equal to 5 or 10 computational grid 
spacings. The posts are tilted at an angle (p to the positive a;-axis, and extend to a height h — Aw above the surface. We start 



considering a rectangular array of posts with lattice constant a^: 
shown in Fig. [2] 



4.6w in the x direction and Uy = Aw in the y direction as 



Placing a sufficiently large drop on the surface to allow significant spreading through the posts is computationally expensive. 
Therefore we define a reservoir of fluid in the centre of the post aiTay typically extending across ~ lOw and reaching to the top 
of the posts. If the fluid density inside the reservoir decreases below the equilibrium density of the liquid phase, new mass is 
slowly added to feed the imbibition. The contact angle of the liquid with both the posts and the substrate is 9 which we vary in 
the range 30° to 70°. This range is representative of different hydrophilic material, for example the polymers considered in|[Tl 
orL16J. 

We first consider a quasi-2D geometry which allows us to concentrate on the directional spreading in the x-direction. We 
choose a simulation box of length AOw in the a:-direction. Aw in the y-direction and Qw along z, with periodic boundary con- 
ditions along both x and y. The reservoir spans the simulation box in y corresponding to simulating a cylindrical drop with 
interfaces which lie, on average, parallel to the ?/-direction. The average of any fluid motion is along x. We then present results 
for a full 3D geometry, using a simulation box with typical dimensions AOw x 32w x Qw, with periodic boundary conditions 
along X and y. A schematic comparison of quasi-2D and 3D spreading is given in Fig.|2^. 



4 

iiKT>f'^ i/fTN' iifn^' 

(a) (b) (c) 




FIG. 3: Advancing front in the unidirectional spreading regime for contact angle 6 = 45° and post angle = 60°. This is a quasi-2D geometry 
with periodic boundary conditions along y. The snapshots (a)-(e) correspond to 0,1,2,4,5 x 10^ time steps. 

III. RESULTS 

Fig.|3]shows snapshots of the imbibition process as a function of time for the quasi-2D geometry and a contact angle 6 = 45°. 
The interface advances in the positive a; -direction, but not in the negative x-direction, because of pinning on the posts. As pointed 
out in Chu et al. |[T] a good understanding of why this occurs follows from assuming that the interface is pinned at the top corner 
of the posts and ignoring any interface curvature along y. We illustrate this situation in Fig.|2]3. The bottom of the interface will 
advance along the surface until it reaches the equilibrium contact angle 6. If this enables the interface to reach the next post it 
will wet this post and move forwards, if not, it will remained pinned because any forward motion will increase the free energy. 
Because of the two dimensional nature of the model, it is easy to work out the threshold Young angles for spreading in the two 
directions as a function of the post geometry. Defining these as 6*+ and 9~ for spreading along +x and —a; respectively, gives 

m 

9,r,±x = tan-1 ( f ) . (7) 

\ax - w T H cot [(I)) J 

Thus there are three regimes: the interface can remain pinned in both directions, advance just along +x, or move forward in 
both directions. The different regimes predicted by Eq. (|7]l are indicated in Fig.|4]as a function of the Young angle 6 and the 
post inclination <j>. Note that, for (p = 90'', 9^ = 9^ as expected. The value of also affects the transition between different 
spreading regimes, as is apparent from Eq.(7): closer posts will ease the spreading, while posts that are further apart will make 
it more difficult. Our choice of — 4.6u' allows us to observe the different spreading regimes over the range of contact angles 
9 and slanting angles (p we consider Fig.|4]also shows the results of simulations for the quasi- 2D geometry. 

All three regimes are reproduced in the simulations. The analytic model gives a good account of the boundaries between 
them, but spreading in both directions is slightly more difficult than predicted by the 2D theory. The simulations allow us to 
identify this as being due to the details of the interface pinning on the posts. For an inteiface advancing along the positive 
X-axis, where the post points towards the direction of travel, we observe three different possible pinning mechanisms, labelled 
\\, B and X in Fig.js] If the post is almost vertical, the leading interface is disconnected and is pinned to the vertical sides of the 
post (I I label). For a more pronounced post tilt the interface remains disconnected, but does not reach the top of the final post (B 
label). For large tilt the leading interface is connected and has reached the equilibrium contact angle on the substrate. However, 
in contrast to two dimensions, the final post is only partially wet. (X label). The situations, B and X, where the interface has 
only reached the top of the penultimate post is only observed for slanted posts. It occurs because the interface can take the 
correct contact angle on the final post without a large penalty in curvature energy. The interface configuration resembles that in 
the partially suspended state identified in Kusumaatmaja and Yeomans H. 
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FIG. 4: Plot of the different wetting regimes in the quasi-2D geometry, as a function of wetting angle 6 and post angle <?!). The other geometric 
parameters are: = 4.6w, h = /Iw. The legend is: blue/circles: bidirectional spreading, red/squares: unidirectional spreading, green: 
no spreading. The different green symbols represent different pinning modes in the forward direction, as described in Fig. |5] Selected 
configurations (black circles) are shown in Fig. [5] The background indicates wetting regimes from a theory assuming a 2D geometry, see 
Eq. m 

We considered a small spacing between posts to facilitate spreading over a wide range of wetting angles. It is worth noticing 
however that, if the spacing between posts were increased in the spreading direction, configurations analogous to the 2D theory 
would likely appear, with pinning on the final rather than penultimate row of posts. 

In the hard direction for spreading, —x, the interface is pinned at the edges of the final line of posts, adjusting to their slope, 
as shown in Fig.|5] label *. Bidirectional spreading only occurs for very low contact angles 9 < 30° or posts close to vertical 
(h > 70°. 



These simulations correspond to quasi-static spreading, with the fluid reservoir replenished very slowly. Borderline configu- 
rations between different wetting modes (such as (p — 70°, 9 = 60°) are very sensitive to exact details of the position and filling 
speed of the reservoir. This is expected because the free energy barriers and capillary forces driving the flow are very small. A 
comparison between two resolutions used shows that, as expected, spreading is slightly more difficult for a narrower interface. 

We also note that, if the rate at which fluid is added to the reservoir is increased, the resulting inertia aUows the fluid to de-pin 
from the top of the posts, forming a spherical cap. 

We next report a full three-dimensional simulation which allows the fluid to spread along both x and y. The reservoir is 
defined as a circular region of radius lOw in the centre of the domain, and the contact angle is 6* = 45°. In Fig.|6]we plot contour 
lines showing the spreading of the drop base for subsequent time steps. The first plot is for vertical posts; as expected it reflects 
the symmetry of the lattice. In the second plot the lattice spacing is the same, but the posts are now slanted with = 60°. In the 
slanting direction the behaviour is consistent with the quasi-2D model, with the liquid spreading only in the positive x direction. 
The spreading in the y direction is comparable to the vertical post case. The overall dynamics closely resembles the imbibition 
observed in experiments by Chu et aI.|[T]. 
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FIG. 5: Different pinning configurations appearing in tlie quasi 2D simulations. The first three correspond to the easy spreading direction, 
X, while the last shows the typical pinning in the hard direction, — x. They can be identified with the circled points in Fig.|4]by matching 
symbols. Cross sections taken through the centres of the posts are shown on the right as full blue lines, (a) If the post is almost vertical, the 
leading interface is disconnected and pinned to the vertical sides of the post (| | label), (b) For more pronounced post tilt the interface remains 
disconnected, but does not reach the top of the final post {B label), (c) For large tilt the leading interface is connected and has reached the 
equilibrium contact angle on the substrate. However, in contrast to two dimensions, the final post is only partially wet. {X label), (d) The 
interface in the negative x-direction is disconnected, and pinned by the sides of the final post (* label). 
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FIG. 6: Three dimensional spreading for 6 — 45° and (a) </> — 90° , Ux — ay = i.Gw. (b) <j> = 60° , Ux = ay = 4.6ty. (c) = 60° , 
ax — 3.6™, ay — 4.6w. (d) (j> — 60° , a^; = 4.6w, ay = 2.53-!/;. Contours showing the interface position at the base of the drop are taken 
every 4 • 10^ time steps. 



In Fig.|6}; the posts are closer in the x-direction. One can see that the asymmetry of the spreading becomes more pronounced, 
with the fluid spreading easily from one row to the next along x, while only spreading very slowly along y. This occurs because 
the thermodynamic driving force for spreading is much stronger in the x direction. For a contact angle of 45° there is slow 
spreading along y, however for a larger contact angle the fluid will remained pinned along y (the crossover can be estimated 
from Eq|7] taking (f) = 90° in Figure 4), and/or x (as seen in the quasi-2D simulations reported in Fig.4). 

Eventually, in Fig. |6}l, we simulate densely packed posts in the y direction, while keeping the same x spacing as in Fig. 
The behavior is now quite different, with the spreading happening first in the transverse direction, and only subsequently along 
the slanting direction. The unidirectionality is however maintained. 

These results indicate that the final shape of a spreading drop can be tuned in detail by varying the lattice geometry and tilt 
angle of the posts, while retaining the relevant property of unidirectionality. 

IV. CONCLUSIONS AND OUTLOOK 

We have applied a lattice Boltzmann algorithm for two phase flow to model the spreading of liquid drops on a surface 
patterned by a lattice of slanted micro-posts. Gibb's pinning of the interface on the sides or top of the posts led to unidirectional 
spreading over a wide range of fluid-substrate contact angles and inclination angles of the posts. Regimes for spreading in no, 
one or two directions were identified, and shown to agree well with a two-dimensional theory proposed by Chu et al| 1 1. A more 
detailed numerical analysis of the contact line configurations enabled us to understand deviations from the two dimensional 
model, and to identify the configurations of the pinned interfaces. 

The final drop shape depends on spacing of the post lattice, the contact angle, and the geometry and inclination of the posts. 
Our simulations coiTespond to slow spreading, but inertial terms will also alter the final drop configuration. Thus there are many, 
varied possibilities to use slanted posts to control drop shapes or the direction of a flowing stream of fluid. Contact angles can be 
varied in situ by electrowetting, and it would be of interest to design substrates with addressable posts where the contact angle 
of each of the posts could be varied independently to allow steering of microfluidic flows. 
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